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Abstract 


of the Jacchia-Roberts (JR), Mass- 

the models was performed to assess the modeling of the total amincnhori^ ^ D J M ) ®* mos P henc density models Evaluation of 
spacecraft and selecting time periods of study representative of all oortions of th^iT h ' S atudy was made generic by using six 
was measured for multiple spacecraft, representing a selection n/^rhi* * ** 11 " yaar solar c V c,e Performance of the models 
altitudes from 400 kilometers to 900 kilometers The orbit □eometh.wt near ' 6quatorial to P° lar inclinations and 

the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) P ^ Earth ' orbitin 9 spacecraft supported by 

employed to minimi^ttit effecte of^mrStng erws" ^TStwt gStentiSmS^ Determination S ' ,sfem ( GT DS) were 

Model-T3 (GEM-T3), was employed to m"nim£ geowteiS S m ^ aVa,la “ e dur,n 9 the anal ^ s . the Goddard Earth 
techniques identified for TOPEX/Poseidon orbit detfrmk^tinn effects on the drag estimation. Improved-accuracy 

atmospheric density riwdd u*^ ^re r |ati U veirhc^^ ne t iC | 3C n Vity near the solar minimum, the choice of 

solar maximum, the differences between the JR DTM and MSIS 86 h 6 " f unn 9 typical solar flux conditions near the 

solar activity, those in which the daily anTei^ay lan Jaf tx tT L° become a PP are nt. Time periods of extreme 

between the models During periods of high geomagnetic activity the standard JR ^ rapi ^’ result ln significant differences 


1.0 Introduction 

Orbit determination for spacecraft whose perigee heights are less than 2000 ,v x 

atmospheric density model because aimncnkprip a er . °°° *" on ieters (km) requires a comprehensive 

these dlimde. 7™ ■? , * 5““ .'"" si 8 nif,c “ l P“‘"' b "™ forces on spacecraft 

atmospheric density models! <GTDS> pr0Vid ' s “»."«» "Mt two 

the Jacchja <*> ^ ^ «» updated to rfpec, 

(Reference 4). At the current rime thfre “L ■ a '' " mod,fed Harris-Prie.ler (HP) model 

Flight Dynamics by the Goddard Space Flight Center (GSFC, 

en^V) I ar^ S the^a^SpectromcteMnco^renF&aner^86 (MSlS-86) ,R Refcren' enl 6) ra, t llre T" “>™> 

have been constructed based on data nna.aiiahic to the JR and HP modi'and 
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vtuying conditions of solar, geom.gn.tic, - 

atmospheric models to determine their potent, al idea tn ' “PP»"' » ^ desin>bk „ les , ttles e models in GTDS 

Dynamics Facility (FDF), particularly in fu rt mi P ' ( di g of tbe GTDS solution fit parameters, 

whme ». evaluation of model perform™. Cation fctom, and definitive ephemeris 
^^“r^ruLlar L geomagnetic conditions, orbit geometries, and spacecraft - «-<* 

ballistic coefficients. 


1.1 Overview . , 

GTDS i, used fo, current opera, ionri orbit ^ ^ ^ 

squares algorithm which estimates the se o or i between observed and calculated values of selected 

parameters that minimizes the sum o « S <1 ?) Qf interest f or this study is the ability of GTDS to 

tracking data measurements over a solution arc (Reference ). 

estimate a drag scaling parameter p, , defined by 


Fn = -~C D pV\V\A{\ + Pi) 
2 


where 

P 

P l 

C D 
A 
V 


density of atmosphere surrounding the spacecraft 
drag scaling parameter, here assumed to be a constant 
: coefficient of drag 
: cross-sectional area of spacecraft 
velocity of spacecraft relative to local atmosphere 


Y — “ y vivviy ~ , 

Assuming th.t the solm flu* and geomagnetic index <GMI) i 

calibrated, then the estimated p\ values shou e near zer ^ assumed to be the most accurate model. A 

sity model that yields the smallest average Pi v ue ^wou ^ give n set of solutions should result in 

model that accurately describes the density magni comparison values would represent poorer model 

minimal definitive ephemens overlap compansons.^ Hg P ' ^ variations over given set of solu- 


1 2 Summary of the Models 

Atmospheric models ate fotmula.ed using "T 

The ora o, semiempirical models began with the ^ 

prime data were derived from atmospheric drag on sateU.tes by use of empirical 

model derived by integration of the d. fflwcw of atmospheric composition data from the mass 
formulas. The MSIS senes of models began with the analyse or a i p ^ q{ ^ ^ wjth data 

spectrometer onboard the Orbital Geophysica serva MSIS-77 model (Reference 10) eventually 

derived from a ground-based radar temperature profile with 

= 52 = — — 
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Table 1. Features of the JR, MSIS-86, and DTM Models 


Model Features 

Jacchia -Roberts 
(Updated to Jacchia 
1971 Standard) 

MSIS-86 

DTM 

Variations with solar flux 

Daily and 81 -day 
centered average 

Daily and 81 -day 
centered average 


Geomagnetic heating delay 

6.7 hours 

0-59 hours, 
with local variations 


Density for each constituent gas 

Only for total 

Yes 

Yes 

Local variations in density 

1 No 

L- — 

Yes, by constituents 

Yes, by constituents 

Diurnal variations 

Yes 

Yes, by constituents 

Yes, by constituents 

Semiannual variation 

Yes 


Yes. bv constituents 

Seasonal latitudinal variations of lower thermosphere 
boundary 




Seasonal latitudinal variations of helium 

Yes 



Hydrogen effects (important above 1200 km) ~ 

No 

Yes 

Yes 


altihiHp^ 3 h ^ CSt ^ ata covera 8® is the 150-km to 600-km range. Jacchia does have some data from higher 

several a , 7’ k * data ° Utside this are ve T sparse, and density extrapolations based on 

several assumptions become increasingly inaccurate for all models as altitudes vary from this region. 

1.3 Direct Comparison of the Models 

s^rTnH riC denSkie , S PrediCted ^ thC JR ’ HP ’ “ d MSIS ‘ 86 mode,s at various a,titudes for common input values of 
r and geomagnetic activity were compared for June 22, 1992. The solar activity is characterized by the 

10.7-centrmeter solar flux, des.gnated in this paper as F )0 7 (with implied units of 10 ^ watts per meter^ per hertz) 

"T? “ by the K, and A„ geomagnetic indices (unitless, where f values a^e 

die locaf T 3 f " r 0 SC ValUCS arC ex P ressed ™ a linear scale). Analysis (Reference 1 1) revealed that 

the local atmospheric density vanat.ons between the MSIS-86 and JR models ranged between -40 percent and + 80 

100> = 2 ^ db — -40 Percent and' +40 percent for 

i 0 ? - 2). These differences represent many subtle differences among the models, the most significant 

local variation being in the behavior of the diurnal bulge. As altitudes increase the diurnal bulgTlv^ntt 

fhitud ^ f° Uth X! l Pri r nly PeIated ‘° ‘ he seasona l _ lat itudinal helium effect, which genejy dominates at 
altitudes above 500 km, where helium flows toward the winter pole. at 

DercenT'duri betWeen f ,ic global averages of densities produced by the JR and MSIS-86 models are less than 15 
percent during low solar activity time periods (F 10 . 7 = 100, Kp = 2). The differences are most pronounced for the 

(F.oT = T 2wT= ifrS'/T values just under 25 percent were observed. At higher solar flux values 

fF - ->nn K* ^ 2 tltudes above 400 km - MSIS-86 predicts smaller densities than JR. Increasing GMI 

(ho.7 = 200, Kp = 5) results in a similar effect. 8 1 

Average densit y plots as a function of time (for altitudes of 300 km, 700 km, and 1300 km) show marked differences 

between the models, particularly in response to GMI fluctuations. The JR model employs a single 1C value 6 7 hours 

pnor to he current time, whereas the MSIS-86 model employs 2. 3-hour A p values spanning^ 7ours prior to The 

red time. Comparison of the time-dependent densities reveal that MSIS-86 and JR are similaMn density 

— * ,OWe ' “ d ™ ddle altitudes for moderate GMI activity but tend to have the strongest reactions to 

dlt, H aC : V y "V m ' dd e ‘° h,gh akltudes - At h '8 h altitudes, JR densities are always greateT than MSIS-86 

ld or e M9 U ,V«°, a , SP i f "h" 6 thC tCndenCy " f ° r ^ exhibi * a ra P' d andLge peak d'shy spike 

and for MSIS-86 to display a broad-based spike. The timing of the MSIS-86 density peak precedes the JR p^lak bv 

^hT M^lT 87 model o" ^ res P onse to solar geomagnetic activity was greater than that of 

rn„HV S 6 .r ' Overall, the largest difference in the models was in their reaction to high GMI activity 

oroduc^dh 4 h d U f? ' r reaC ! 10nS ‘° the SO,ar flux aIso differcd The >arge differences observed in the densities 

H I ^ T !* a significant effect on areas such as mission planning, where the density is 

important in determining orbit decay rates. e uensuy is 
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Orbit Determination Parameter or 
Option 

Estimated parameters 


Table 4. Parameters and Options Used in the GTDS Solutions 

Parameter or T -- — — 

User Spacecraft TDI 


Integration type 


Coordinate system of integration 


Integration step size (seconds) 


Tracking measurements 


Data span 


Data rate 


Editing criterion 


Measurement weight sigmas 


Satellite area model (all constant) 


Geopotential model 


Atmospheric density model 


Solar and lunar ephemerides 


Coefficient of drag (Cp) 


User-spacecraft antenna offset 


Tropospheric refraction correction 


Ionospheric refraction correction 
Ground-to-spacecraft 
Spacecraft-to-spacecraft 


Polar motion correction 


Solid Earth tides 


Orbital state, p, , and station measurement 
biases (USO bias and drift, COBE only) 

Cowell 1 2th order 


Mean of J2000.0 


60 seconds 


TDRSS two-way Doppler (TD2S) 
TDRSS two-way range (TR2S) 
TDRSS one-way Doppler (TD1 S) 
Ground S-band range rate (URDF) 


2 days (4 for COBE) 


1 per 1 0 seconds 


3o 

Central angle to local horizon 


TD2S: 0.25 hertz 
TR2S: 30 meters 
TD1S: 0.13 hertz 
URD F 1 0 centimeters/second 


Orbital state, coefficient of solar radiation 
pressure (C R ), BRTS range bias 


Cowell 1 2th order 


Mean of J2000.0 


600 seconds 


BRTS two-way range 


See text 


1 per 1 0 seconds 


3a 


10 meters 



COBE: 17.8 meters 2 

ERBS. 4.7 meters 2 

HST: 74.0 meters 2 

LA4: 12.3 meters 2 

LA5 12.7 meters 2 

SMM: 17.5 meters 2 

40 meters 2 

COBE. 2155.00 kilograms 

ERBS: 21 16 00 kilograms 

HST: 1 1 328.00 kilograms 

LA4 1900.32 kilograms 

LAS: 1913.25 kilograms 

SMM: 2315.59 kilograms 

TDRS-4 - 1900 kilograms 
TDRS-3 - 1990 kilograms 
TDRS-1 - 1730 kilograms 

50 x 50 GEM-T3 

20 * 20 GEM-T3 


JR, MSIS-86, DTM 


_DE200 


2.2 (2 3 for COBE) 


Constant radial 
Yes 

Yes 

No (central angle edit instead) 


N/A 

DE 200 


N/A 

No 

Yes 
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r 100 



Day of January 1 987 

Figure 1. Solar Flux and GMI for Period A 


Table 5. Summary of Orbit Determination Solution Results for Period A 


Spacecraft 


ERBS 


Landsat-5 


SMM 


Density Model 


JR 

MSIS-86 

DTM 


JR 

MSIS-86 

DTM 


JR 

MSIS-86 

DTM 


PI 


-0 01 ± 0.22 
-0.07 ± 0.25 
-0.24 ±0.18 


-0.01 ±0.31 
-0.01 ±0.31 
-0.16 ±0.26 


-0.27 ± 0.04 
-0 14 ±0.03 
-0.37 ± 0.02 


WRMS 


0.21 ±0.04 
0.21 ± 0.04 
0.21 ± 0.04 


0.21 ± 0.05 
0.21 ±0.05 
0.21 ±0.05 


0.26 ± 0.03 
0.24 ±0.02 
0.25 ± 0.02 


Overlap MPD (meters) 



11/20 11/21 11/22 11/23 11/24 11/25 11/26 

Day of November 1 990 


Figure 2. Solar Flux and GMI for Period B 
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23.5 ± 6.6 
23 5 ± 6 6 
23.5 ± 6.6 


27.7 ±10.9 
27.3 ± 10.3 
27.3 ±105 


26.7 ± 7.8 
23.9 ± 9.2 
24.6 ± 9.3 
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EpochDate 

Figure 6. ERBS Solution WRMS Results for Period C 



Table 7. Summary of Orbit Determination Solution Results for Period C 


Spacecraft 

Density Model 

P\ 

WRMS 

Overlap MPD (meters) 

COBE 

JR 

MSIS-86 

DTM 

DTM-C 

-0 21 ±0.04 
-0 01 + 0 06 
-0 8210 34 
-0 1410.15 

0 22 1 0.02 
0 231001 
0.36 1 0.07 
0 32 ± 0 04 

30 4 1 12 4 
41 8 1 17 2 
103 4134 1 
82 4 + 21 8 

ERBS 

JR 

MSIS-86 

DTM 

DTM-C 

-0 0810 05 
-0 04 1 0 08 
-0 50 + 0 38 
-0.0710 18 

0.20 1 0 05 
0 22 ± 0 09 
0 3310 16 
0.29±0 13 

22 4 ± 9 3 
33 3 ± 17 7 
77.6 ±37.9 
61 2 ±316 

HST 

JR 

MSIS-86 

DTM 

DTM-C 

-0 24 ±0 05 
-0 13 ±0 05 
-0 27 ± 0.26 
-0.20 ±012 

0.28 ±0 10 
0 28 ±0 09 
0.55 + 0.31 
0.45 ±0 24 

39 6 +24 0 ~~ 

50 8 ±36 6 
115 51845 
85 1 ± 70 3 

Landsat-4 

JR 

MSIS-86 

DTM 

DTM-C 

-0 23 1 0 04 
-0 00 + 0 10 
-0 6010 42 
-0 03 + 0 20 

0 17 + 009 
0 19 + 0 10 
0.2610.15 
0 24 + 0 12 

21.5 + 79 
37 9 ±23 9 
87 4 ±45 4 
69 3 ± 32 2 
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OvenUI. 'he m and MSIS-86 models were comparable In performance, while both versions of 'he DTM model 
milled the worst WRMS «d overlap MPD values, on average All the models had some difficulty wiiVhe CM 

S 3 modllsT r-T n 5 ‘ “ WOrSI Were ' m - MS,S ' 86 ' D ™‘ ■* -ltd implemented DTM 

dS.fLn , ^ boot DTM versions produced overlap 

FguITs) m " r “ ,Uir " n ' nl) for HST «e» CMI activity on February I (see 


*♦- DTM -A MSISaM-IW- DTM-C 



Figure 8. HST Solution Overlap MDP Results for Period C 


The p } values resulting from the originally implemented DTM model based on the 81 -day endpoint average solar flux 
are significantly higher then the JR and MSIS-86 p x values, while the modified DTM model using the 81 -day 

“ n,er . J c 7™ g ! S ° Iar , flUX reSUUed in aVCrage P[ com P arab ' e “> ^ and MSIS-86. Daily trends for the solved-for 

JLT V fj 2 ^ S afe f nearly C ° nStant - 1116 D ™ Pi values slart 0, peak at approximately 0.7 for HST 

(higher for the other spacecraft), and then return to the -0.2 to 0.2 range. Tl,e modified DTM model showed a 

smnlar peak; however, the relative height of the peak from the base was onjy half that seen in the implemented 
I. The atmospheric density modeled by the DTM model was consistently low during the peak in the daily solar 

^rtion of tb h ' nd,Cated fl by thC A , tre " ds in both of ,he D ™ cases This is consistent with a situation in which a 
portion of the daily solar flux is applied as a daily difference from the mean, as opposed to being correctly applied as 

part of the mean value^ The mean solar flux has a greater effect on the resulting density than the daily contrition 

solar n fi. h, H “ 8, * d V endpoint-averaged F 10 7 value will result in a significant portion of the 

solar flux during the peak activity being applied as the daily solar flux input, resulting in a lower density. The fact 
that the phenomenon is still apparent with the modified DTM indicates that the DTM model may not handle extreme 
solar flux input values as well as the JR or MSIS-86 models. 


3.4 


Period D: Effects of Geomagnetic Activity 


Period D was chosen from a time period in which the daily and average solar flux values were behaving in the 
nominal 27-day period pattern, as in period B, but also included extremely high geomagnetic activity. As shown in 

f l 3Verage S °‘ ar , X , W3S appr0Ximatel >' 200 > while the daily value was approximately 240, near the 
maximum of the current 27-day solar rotation. In general, GMI activity was very high, with K p actually reaching 8 7 

(9.3 ,s the nominal maximum on the logarithmic Kp scale) on one occasion. The behavior of the atmospheric models 
dunng such geomagnetic storms is important because the ability or inability to accurately model GMI activity effects 

THR^Tn J 6 °?m determina ‘ ion and prediction process adversely, as has been observed during the operational and 
TDRSS Onboard Navigation System (TONS) experiment use of the JR model (Reference 13). 

BaS ed on the generally good performance of the DTM model during this period, a second set of results was generated 
for the JR model, in this case, the delay in the geomagnetic activity was modified from the original 6.7 hours to 
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1 hours This change was made based on previous analysis of the performance of the models (Reference II), which 
have®, he longest delay in OMI, and because ,he value « 

assumed value meant to reflect an average time for the geomagnetic heating (R 

model is referred to as the JR-3 model in this paper. 



Figure 9. Solar Flux and GMI for Period D 


As is evident from Table 8, there were significant variations in the p\ values. The p\ values ^suiting ™ 

r“ ^ U ^«ERBsX~Tfl .ri re active <L to, oih.r grinds, due in .he geomagnetic activity^ In this 
case.^he MSIS-86 model appears to produce a consistently lower average density than either JR or DTM or t ose 
spacecraft that are under 700 km. 


Table 8. Summary of Orbit Determination Solution Results for Period D 


Spacecraft 

Density Model 

Pi 

WRMS 

Overlap MPD (meters) 

COBE 

JR 

JR-3 

MSIS-86 

DTM 

-0.28 ±0.06 
-0.28 ± 0.04 
-0.18 ± 0.06 
-0.16 ±006 

0.25 ± 0.05 
0.23 + 0.02 
0.20 ± 0.02 
0.21 ± 0.01 

49.1 ± 19 6 
30 9 ± 11.5 
28 0± 6.6 

19.1 ± 59 

ERBS 

JR 

JR-3 

MSIS-86 

DTM 

-021 ±0.10 
-0.20 ± 0 06 
-0 06 ± 0.09 
-0.27 ± 0.06 

0.22 ± 0.08 
0 20 ± 0.04 
0 23 ± 0.08 
0.22 ±0 06 

58.7 ± 56 4 
33.0 ± 27.2 
53.3 ±43.2 
35.8131.5 

HST 

JR 

JR-3 

MSIS-86 

DTM 

-0.37 ± 0 09 
-0.38 ±0.07 
-0.27 ± 0.09 
-0.52 ±0 05 

0.55 ±019 
0 45 + 028 
0.50 ± 0 25 
0.35 ± 0 24 

62 9 ± 47.0 
49.8 ±26 3 
77 4 ±60 5 
51 2 ± 26.8 

Landsat-4 

JR 

JR-3 

MSIS-86 

DTM 

-0.29 ±0 10 
-0.38 ± 0 07 
-0.16 ±0.15 
-0.29 ± 0 06 

018 ±0.05 
0.16 ±0.04 
0.17 ±0 03 
017 ±0 04 

36 2 ± 13.2 
31.1 ±22.2 

37 5125.7 
33 41 19.4 
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[-B-JB ■♦■DTM -A- MSS«-»-JB.3 ~~| 



Figure 10. ERBS Solution Results for Period D 


The increased WRMS and MPD values during the GMI storm on June 4 through June 7 show that there was 
sensi .vity to the GMI, especially for the models with the longer delay in GMI modeling (JR and MSIS-86) Solution 
WRMS values showed significant variations between the models and were adversely affected by the GMI activity 

GMI D8 Jt r , SU ~' thC COBE WRMS ValUeS from the D ™ model we« little affect* by the chmge i^ 
GMI activity on June 6 and June 9, unlike the JR and MSIS-86 models, both of which showed increased WRMS 
values. The modified JR produced individual solution WRMS values similar in trend to DTM. 

For all models the worst overlaps are seen at the onset and ends of the large storm from June 5 through June 9 The 
MPDsare^ce^ “it MSIS ' 86 models P roduced over >a P s of approximately 175 meters in the 2a to 3a fange The JR- 3 

°7 th : Standa 1 * modeI 1116 Landsa ‘- 4 «"**P values are significantly lower than 
. f 3 Ues for l ^ ls P erlod due to the absence of the overlap for the June 4 and June 5 solutions 

™ ZTT ms r icd shows - improv '' m ''' whe " ' he '■» i” 

MSIS 86 ’ WeVef ’ IS n0t f ‘ hat jt U the best P erformer when compared with DTM and 

3.5 Long-Term Changes in Density Model Performance 

The long term behavior of the estimated p, values is of interest because it indicates long-term variations in the 
modeling of the atmospheric density. To do this accurately, it is necessary to consider fhe Z va ue foVthose 

COBE 10 lhiS *“* ,he VZ 

™ c S d C0BE F,gUre illustrates the average />, values (Rhol in Figure II) for each model 

used in the ERBS orbit solutions. (JR -3 applies to period D only, while DTM-C applies to period C onlvl The n 

t V he“S p is°f ^ ^ MS ,' S n 6 m ° ddS Change Signif,Cantly de P end ' n g on the study period. The total ra^e of 

varies abo buTtheTotal ' *° reprCSentlng up *° 67 P ercent of the actual atmospheric density. DTM 
fl|1 b * t0t f’ range . ,s somew hat smaller assuming that the DTM using the 81-day centered average solar 
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Period 


Figure 11. ERBS Long-Term Solution p-, Changes 


4.0 Summary and Conclusions 

Performance of the three atmospheric density models was measured for multipj. spacecm*. 'XTd'tclL“onsof 

* - 

GSFC FDD. 

=3 E=iErsxrrs== 

cases the WRMS values and the overlap MPDs result in similar conclusions, 
models. Generalfy. the IR model performed the best, while DTM perfomed the worst. 

original paper for optimal model performance. However, the improvements ,n the DTM performance resulting rom 
this change were not enough for the model to outperform the JR model. 

th^smrn^d^R'TO^h'^^c^a^a'd^^our^U^for'gm^^nct^^ffccl^was^tper^nimd 8 ^ DT^.|wh^r'^a^ 

effects is particularly critical. 
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{££ "1 ? C rCnCCS m ‘ he ° rbit determinatio " performances of these models is negligible 

ornfm»l n f un s ett ^ ge°magnet ,c activity, the JR model currently implemented in GTDS did no. provide 

optimal performance. With the exception of COBE, the DTM model appeared to handle GMI activity best for the 
spacecraft used dunng that period. Modification of the JR model geomagnetic activity modeling to reflect a 3-hoSr 
delay instead of the default 6.7-hour delay produced results that were similar to or marginally better th^ the DTM 

is “ in — 
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